Josephson current through a quantum dot with spin-orbit coupling 
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The equilibrium Josephson current through a nanoscale multi-level quantum dot with Rashba 
or Dresselhaus spin-orbit coupling a has been computed. The critical current can be drastically 
modified already by moderate a. In the presence of a Zeeman field, Datta-Das-like oscillatory 
dependencies on a are predicted. 

PACS numbers: 74.78.Na, 74.50.+r, 71.70.Ej 

I. INTRODUCTION 

The Josephson effect provides a fundamental signature of phase coherent transport through mesoscopic samples 
with time reversal symmetry [l|. Here, we reconsider the theory of the Josephson effect through a quantum dot (QD) 
in a two-dimensional electron gas (2DEG), taking into account Rashba and/or Dresselhaus spin orbit (SO) couplings 
0, Despite the large recent interest concerning SO effects in QDs, mainly caused by quantum information Q 
and spintronics applications [f| , the question of how the Josephson current is modified by SO couplings has barely 
been addressed. Naively, since SO couplings do not break time reversal symmetry one may expect that they will not 
affect the Josephson effect at all. This expectation seems basically confirmed by the existing theoretical studies: In 
Ref. [H , the Josephson current through a perfectly contacted 2DEG with Rashba coupling is predicted to only show 
SO-related modifications in the simultaneous presence of a Zeeman field. In Ref. Q, the same conclusion has been 
reached for a single-channel conductor with arbitrary contact resistance. Moreover, in Ref. j|, by a generalization of 
earlier scattering approaches [§, [nj, the case of a wide 2DEG with arbitrary contact resistance was studied, where 
the semiclassically averaged Josephson current is again found to be insensitive to Rashba SO couplings. However, 
the presence of a Josephson current through the mesoscopic system necessarily breaks time reversal invariance and 
therefore we reexamine this expectation in the present paper. 

As 2DEG devices based on InAs-related materials are known to exhibit strong gate-tunable [ld[ Rashba SO 
couplings, and possibly Dresselhaus SO couplings, we study this problem here from the perspective of multi-level 
QD physics. Sup ercurrents through related devices have been probed by several experiments; for recent work, 
see Refs. [lU, Il3l . ITil . flH [l6j . It is thus not only of academic interest to quantitatively examine the effects of 
Rashba/Dressclhaus SO couplings on the equilibrium Josephson current. Moreover, very recently, gate-tunable su- 
percurrents through thin InAs nanowires have been reported (l7l [l8| , revealing complex current-phase relations such 
as 7r-junction behavior. Although we study a 2DEG geometry, our results are also relevant for such nanowires: the 
transport channels reside in a surface charge layer, and Rashba terms due to narrow-gap and strong confinement fields 
dominate over all other SO couplings. In those experiments, Josephson currents through few-level dots have already 
been achieved. 

In the present work, we arrive at the surprising conclusion that SO couplings have huge effects on the Josephson 
current in nanoscale multi-level dots. Depending on the parameter regime, the critical current can be largely sup- 
pressed or enhanced. In a single-level dot, SO effects on the Josephson current vanish unless there is also a Zeeman 
field, but in multi-level dots no such restriction applies. The predicted strong dependence on SO couplings should be 
observable in state-of-the-art experiments. We also show that oscillations in the critical current I c as a function of 
the distance L between the lead contacts appear under suitable conditions, resembling Datta-Das § spin precession 
effects. 

The structure of the remainder of this paper is as follows. In Sec. [Hj the model is described. The Josephson current 
is calculated in Sec. IIII1 and results are presented in Sec. IIVI for the simplest cases of one and two spin-degenerate 
levels. We conclude in Sec. [V] Some details concerning the derivation of the Josephson current have been delegated 
to an Appendix. 
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II. MODEL 



We study a QD formed by a confinement potential V(r) within a 2DEG, r = (x, y). The total Hamiltonian reads: 

H = H D + H T + H L +H R (1) 

and using the QD fcrmion creation operator d£.(r) for spin a = db, the isolated QD is described by (we put h = k R = 1, 
and spin summations are often left implicit) 

H D = jdr dt(r) ( { ~ lV + *l~ 2a2 +*.9+v) d(r), (2) 

where d = (df, c?^ ) T ' , m is the effective mass, b = (6 X , b y , b z ) is a constant external Zeeman field (including gyromagnetic 
and Bohr magneton factors), and a = (a x , a y , a z ) with standard Pauli matrices. Orbital magnetic fields can also be 
taken into account in our formalism but give no qualitative changes. In Eq. @, the x(y) component of the operator 
a = (a x , a y ) acts in spin space, 

a = a ( sin 9 a x — cos 9 a v , cos 9 a x — sin 9 a y J , (3) 
and contains the combined effect of Rashba (a R ) and linear Dresselhaus (od) couplings via 



a=ya R + a D , sm9 — ao/a. (4) 

These two are generally the most important SO couplings in QDs based on 2DEG geometries. The superconducting 
banks are described as 3D s-wave BCS models. The Hamiltonians of the left and right superconducting electrodes 
have the standard BCS form, 

H j=L/R = ^k^^ + Y,{ AeTlH2 *ki*l(-vi+ h - c )> ( 5 ) 

with £k = k 2 /(2m) — fj, and electron creation operators ^ L / R ^ a , and carry the relative phase 4>. For simplicity, 
we assume the same BCS gap A on both banks; in experiments, both contacts are typically created by the same 
lithographic process. Finally, Ht describes tunneling between the leads and the QD, 

H T = t L/R,n 1 ^[/ RM U L / R : d n + h.c, (6) 

L/RM,n 

where for now U^/ R = 1 (but see below), and the index n labels the QD eigenstates for a = b = 0, 

^V 2 + V(r)) X «(r)=e„Xn(r), (7) 

with d a (r) = ^ n Xn{ v )d n a- To simplify our expressions, we make the inessential assumption of k- and spin- 
independent real- valued tunnel amplitudes t L j Rn between the L/R banks and the QD state \ n . Although the 
connection between the superconducting leads and the dot appears in the form of a tunneling Hamiltonian, we treat 
these terms below to all orders by integration over the lead degrees of freedom. This provides justification for ig- 
noring Coulomb interactions on the QD. Over the last decade, several works [l9[ have shown that in the limit of 
high transparency, i.e. for tunneling rates large compared to the dot charging energy, charge quantization (Coulomb 
blockade) effects are strongly suppressed due to large charge fluctuations. Note also that an observable Joscphson 
current usually requires high-transparency contacts. 

Without loss of generality, we now choose a coordinate system where the tunnel contacts are located at (x, y) = 
(=FL/2, 0). It is then convenient to eliminate the a^-term by a unitary transformation on the dot fermions [201 ] . 

d(r) e- lxa *d{r). (8) 

When expanded into the basis ((7|), the transformed Ho is 

H D = ( e n - d n d n + 4 ( A «n' + B nn /) • 5dn>, (9) 
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where the transformed a y is encoded in 



cos#[l - 2cos(26») sin 2 (ax)] 
A w / drx*(r)[fy Xn ,(r)] ( - sin fl[l + 2 cos(20) sin 2 (a;r)] | , 

— cos(26>) sin(2ax) 

and the Zeeman field leads to 

Bw • B= J dr X * n (r) (b • e***-^-***-) Xn ,(r). (11) 

Since the dot electron wave function x«( r ) in Eq. ([7|) can always be chosen real-valued, the Ai^ nn i (-Bj, ran ') with 
i = x,y,z represent antisymmetric (symmetric) Hermitian matrices in the QD level space. Finally, the transformation 
© yields the 2x2 spin matrices 

_/ cos(aL/2) Tsin(aL/2)e- ie \ 
"l/h- ^±Bin( a L/2)e* cos(aL/2) J 

entering Ht in Eq. The oscillatory dependence on the length L between the tunnel contacts appearing for a^O 
is due to a spin precession phase [f|. 

III. JOSEPHSON CURRENT 

The equilibrium Joscphson current at temperature T follows from the free energy —T\nZ via 

I{4>) = ~T8^\xxZ{4>), (13) 

where the superconducting phase difference enters the BCS Hamiltonian ([5]). We adopt a functional-integral repre- 
sentation of the partition function Z , which here requires a slightly non-standard formulation due to the spin-flip 
terms. Given the fact that the Hamiltonian is quadratic in the fermion operators, the calculation should proceed 
in a straightforward manner. In the absence of spin flips, such quadratic Hamiltonians are typically dealt with by 
introducing Nambu spinors. For the present case (see Appendix), the presence of spin flip processes makes it necessary 
to introduce two types of Nambu spinors to describe the dot. In the end, the effective action describing the quantum 
dot coupled to the leads can nevertheless be written in terms of Grassmann fields d n<T (r) and d na (T) for the dot 
fermions (r is imaginary time). For a multilevel dot, we thus form the Nambu spinor for dot level n, 



d_n]{r) 



= Vf£V fa,T A l (w), (14) 



with fcrmionic Matsubara frequencies uj = (21 + 1)ttT (integer I). Similarly, we define the Nambu multispinor 

D{u) = {D 1 {u),D 2 {u),...) T , (15) 
and then integrate out the lead fermions. As a result of this operation, Z = J T>(D, D)e~ s with the action 



(D^^Di-^S^P^y (16) 



The apparent doubling of spinor space reflects our above discussion. This doubling does not create a double-counting 
problem, since only positive Matsubara frequencies lo appear in Eq. (|16[) . After straightforward but lengthy algebra, 
with A± inn i — A x<nn i ± iAy :nn i (and likewise for B), see Eqs. (|10[) and (fTTj) . we find 

a _ f A z g z + B z A^g x + iB-Oy \ 

u ~\ A+a x -iB +( r v -[Gzlf + A z a z - B Z J ' [U) 

The 2x2 structure refers to this doubled space, where each of the four entries represents a matrix in multispinor 
Nambu space, i.e. the Pauli matrices in Eq. (fTT)) act in the Nambu (spin) sector, while A and B operate on the level 
space. Spin-flip transitions caused by SO processes are described by the off-diagonal entries A± and B±. Finally, 
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Guj denotes the QD Green's function after integration over the leads but in the absence of SO couplings. This 
Green's function is derived and specified in the Appendix using the wide-band approximation, see Eq. (|A2[) . with the 
leads characterized by a constant normal-state density of states v (assumed identical in both banks), and with the 
hybridization matrix 

^L/Rjin' = Kv\t L / Rn t L / R . n ,\. (18) 

Note that eventually the Vl/r. factors do not contribute to the action, see Eq. (fl"7|) . and thus the precession phase aL 
can appear only via A and B. We finally obtain from Eq. (|13j) the equilibrium Josephson current-phase relation, 

I(<t>) =~Tj2^lndetS u . (19) 

ui>0 

We stress that this result (given the model) is exact and does not involve approximations. The calculation of the 
Josephson current requires only a simple numerical routine and allows for a quantitative comparison to experimental 
data. 



IV. APPLICATIONS 



In order to get insight into the relevant physics, we here analyze the simplest cases of a single and of two spin- 
degenerate QD levels with symmetric hybridization matrix, r£ jrm < — Tji, nn i = T nn i/2. While our conclusions are 
general, some coefficients below are model-dependent. They are derived for V(r) given as hard-box confinement 
along the x-axis, —L/2 < x < L/2, plus a harmonic transverse confinement of frequency scale lu±. The level index 
n = (n x ,n y ) then contains the respective integer quantum numbers n x > 1 and n y > 0, with eigenenergy (up to an 
additive constant related to a gate voltage) 

e„ = (7r ^ L) + u x (ny + 1/2) (20) 

in Eq. (|7|). Note that Eq. (fT9|) contains not only the contributions of Andrcev bound states but represents the full 
Josephson current. The necessity of computing the full current for finite L has been stressed recently [2lj |. 



A. Single dot level 

Let us first address a single QD level, where we reproduce and extend previous results @, 0] ■ The antisymmetry of 
A„„' implies A = 0, and thus SO couplings do not affect the Josephson current in the absence of a Zeeman field [f|. 
With B = |B| given by 

B 2 = (b x sin 6 - b y cos 9) 2 + F„ x [b\ + {b x cos 6 + b y sin 9) 2 ) , (21) 

where F n is denned as 

sin(aL) 

the determinant entering Eq. (fll)j) is 

, a 99/9 9 r 2 A 2 cos 2 (<i/2) 
detSu = 4 7 2 B 2 + f e 2 + 7 2 + ^ 2 - ^ ' - B 2 ) , 



F n {ot-L) = [ aL{1 _{ aL j to)2) ) , (22) 



where 7^ = ui(l + T/^lu 2 + A 2 ). Here we consider a situation where transport proceeds through a single resonant 
level with longitudinal quantum number n x > 1, i.e. e = — a 2 /2m is close to zero, while the dot stays in the 
transverse ground state n y = 0. Under the condition ui± ^> Se 3> |b| for the level spacing 5e, i.e. for a quasi-lD 
situation, all states besides the n x state close to resonance can be neglected, and one effectively has a single-level dot. 
Clearly, SO effects now can only enter via the effective magnetic field (|21[) . Compared to the Andreev level spectrum 
of the a = b = QD connected to superconducting electrodes, the combination of SO and Zeeman couplings gives 
rise to a splitting of the Andreev levels, giving in total four distinct levels. Each level carries a definite spin with 
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both spin up/down levels above/below the chemical potential. Superficially, the structure of this Andreev spectrum 
is somewhat analogous to the one of Ref. [H}. However, it is actually different because in Ref. [13] the levels refer to 
zero spin, spin up, spin down, and a spin singlet state, respectively. With Eq. ([23]) the Joscphson current appears in 
precisely the form of the mean-field solution of the interacting Anderson dot problem obtained by Vecino et al. [23| . 
This implies that their results can be taken over (with their exchange energy E ex being our B). In particular, the 
phase diagram in the B-e plane exhibits rich behavior with four different phases including 7r-junction behavior, see 
Fig. 3 of Ref. [23j . In Fig. [TJ the Joscphson current is plotted as function of the phase difference. For simplicity, we 
consider e = and gradually increase the effective magnetic field B given in Eq. (|2ip . The transition to a 7r junction 
with negative critical current occurs when B > s/T 2 + e 2 . 

For b x — by — 0, the oscillatory behavior of B due to the dependence of F n on the spin precession phase aL is 
most pronounced, with B = for aL = 27r/ (integer I). This oscillation may then persist in the Joscphson current, 
suggesting the appearance of Datta-Das like oscillations in the critical current. These oscillations are displayed in 
Fig. [H for n x = 1, such oscillations are not yet observable, but they become visible for n x > 1. Note that for 
convenience, the plots in Fig. [2] ignore the a 2 /(2m) shift of e„ x , as is appropriate for small a; alternatively, one could 
adjust a gate voltage on the quantum dot in order to bring the level into resonance. For T ^ A, |b|, the amplitude of 
the current oscillations is of the order of (eA/ft)(|b|/r) 2 . The oscillationperiod in the critical current is roughly set 
by aL — n, similar to the normal-state case of the Datta-Das transistor Moreover, by systematic variation of the 
magnetic field (b) direction within the x-y plane, one could measure the SO angle from the Josephson current, see 
Eq. 



B. Two-level dot 



In order to get SO-related effects on the Josephson current without a magnetic field, b = 0, we study a QD with two 
levels as a minimal model. As a concrete example, consider the dot states given by the first two oscillator eigenstates 
{n.y = 0, 1) in the longitudinal ground state (n x = 1), such that the energy levels are encoded in e 2 ~ £i = w ± an d 
e o = ( € i + e 2)/2. [Note that n y ^ n' y is necessary for A„„/ ^ 0, see Eq. (fTOj) .] Now all matrix elements of B and A 
are zero except for Ai 2 = — A 2 i, with A = \ A 12 \ given by 

2 

A 2 = ( S in 2 (20) + F 1 (aL) cos 2 (20)) . (24) 

2m v ' 

One easily checks that detS^ and therefore the Josephson current now depends on a and exclusively through 
Eq. (24]). 

Figure [3] shows typical numerical results based on Eq. Q19p for the critical current I c as a function of a for two 
choices of the hybridization parameters T nn , with Ti2 = \/riir 22 . Here we consider a pure Rashba coupling (0 = 0), 
and the dot levels are determined by eo = and hiu± = 20 mcV. The shown range for aL can be realized in InAs- 
based devices (Tlj . and thus the supercurrent can be strongly modified by experimentally relevant SO couplings in 
multi-level quantum dots, with pronounced minima or maxima in I c . The apparent cusp-like maximum in Fig. [3]is 
smooth and does not represent singular behavior. 

Note that the results displayed in Fig. [3] are for large T/A, where the Josephson current is predominantly carried 
by the Andreev bound state contribution Ia{4>)- The latter can be analytically evaluated from Eqs. (|17p and (fT9l) for 
(T\i + r 2 2)/A S> 1, and (at T = 0) is expressed in terms of an effective transmission probability To as in Refs [Hill, 

IaW = -7-r I (25) 
2h ^1 -T sin 2 (0/2) 



where we find Tq = 1/(1 + z 2 ) with 



A 2 +u;l/4-(eo-a 2 /2 m y 



(Tn + r 22 ) (e - a 2 /2m) + (r u - r 22 ) w x /2 ' 



The minimum in I c is thus expected when To ~ 0, in accordance with the values seen in Fig. [3] (For the dashed curve 
in Fig. [3J the I c minimum is at a = 0.) On the other hand, the maximum in I c (where I c may exceed the a = 
value) is found for Tq = 1, again explaining the numerical result in Fig. [3J Here the level shift a 2 /(2m) brings one 
dot level into resonance and thereby causes the enhanced supercurrent. This simple reasoning also explains why the 
position of the I c maximum depends only very weakly on the T nn > . 
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The resonance condition, however, depends trigonomctrically on the SO coupling a via A 2 , see Eq. (|24[) . It results 
from the interplay between SO couplings and the multilevel nature of the QD. By the description in terms of Andreev 
bound states, we can also understand why for large dots (long or wide), the effect of the SO coupling is negligible. 
First, for a 'long' dot, taking 6 = 0, A 2 decreases as L~ 2 , and therefore the SO effects on the Josephson current 
vanish as L — ► oo. Second, when the lateral dimension of the dot is large, the confinement frequency ujj_ becomes 
small. Now A 2 contains tu± as a prefactor, and hence is once again suppressed. In both situations, the effect of the 
SO interaction is merely reduced to a shift of energy levels. 

A qualitative argument in favor of such a resonance process can be reached as follows. Schematically, when a 
Cooper pair enters the dot from the left, either its electrons occupy the same transverse level or they can choose 
different 'paths' (different transverse levels), in a manner quite similar to a cross- Andreev scattering process [Hj. The 
SO interaction acts differently in these two levels because the longitudinal momenta of both states are not identical. 
Two electrons entering the QD from, say, the left electrode have initially antiparallcl spins, but their respective spins 
now precess at different rates because of this mismatch in longitudinal momentum. Depending on the value of aL, the 
two electrons which then exit the QD at the right side may, however, be brought back to an antiparallel configuration, 
leading to a resonance in the critical current caused by the presence of SO couplings. On the other hand, if the spins 
do not reach the antiparallcl configuration, the Josephson current will be reduced. 

Note that more than one resonance can be achieved by a careful choice of parameters. This effect is illustrated in 
Fig. 2J For a multi-level dot with small confinement frequency, the multiple resonances associated with all possible 
pairs of paths (pairs of levels) will start to overlap significantly. It is then natural to expect that in the limit of 
very wide quantum dots, these averaging effects will wash out any SO-related structures in the Josephson current, 
consistent with the results of Ref. @ for infinite width. 

V. CONCLUSION 

To conclude, we have studied the Josephson current through a multi-level quantum dot in the presence of Rashba 
and Dresselhaus spin orbit couplings. In the absence of electron-electron interactions on the dot, this problem is 
exactly solvable, a simple consequence of the fact that the Hamiltonian is quadratic in the creation/annihilation 
operators. (Our model consists of a noninteracting quantum dot including spin flip processes and connected to 
BCS superconducting leads.) Nevertheless, the combination of superconducting leads and spin-flip terms renders 
a calculation of the Josephson current a quite technical task, which we chose to address via functional integral 
techniques. We have explicitly shown results for one and two levels, but our general expressions can be applied to 
arbitrary situations. 

For a single dot level, spin-orbit effects cancel out unless a magnetic Zeeman field is included. In this case, we 
predict spin precession (Datta-Das) effects in the Josephson current, i.e. oscillations as a function of the effective 
length of the dot. These oscillations have amplitude of the order of a few tenths of the nominal critical current, 
which should be observable. They result from the interplay between the period of the oscillating effective magnetic 
field (caused by the combined effects of the Zeeman and Rashba interactions) with the wavelength of the longitudinal 
modes in the dot. 

More interestingly, for the case of a double dot, spin precession effects show up even in the absence of an external 
magnetic field. This is a novel effect in the field of Josephson physics for devices subject to spin-orbit coupling. The 
supercurrent can be drastically modified, either containing sharp peaks or being largely suppressed. The experimental 
observation of such peaks could constitute evidence for spin-orbit effects in a superconducting transistor. 

Possible extensions of this work could include Coulomb interaction effects in the dot, which will be important when 
the tunneling rates become comparable to the dot charging energy. While the inclusion of such effects is beyond the 
scope of this paper, the present formulation of the problem with functional integral techniques can be adapted to 
include them in an approximate manner. Nevertheless, we stress again that the high transparency limit considered 
here is in fact quite relevant in the light of recent experiments where quantum dots are embedded in a Josephson 
setting using InAs nanowires [l?], EH ■ 
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APPENDIX A 

In this appendix, we provide some details concerning the derivation of the action in Sec. [H] The partition function 
is here calculated by using the path-integral approach. For this purpose, we rewrite the Hamiltonian in terms 
of Grassmann-Nambu spinors. We introduce two types of Nambu spinors ('spin-up' and 'spin-down') for the QD 
fcrmions, 



while the BCS leads are described by the conventional 'spin-up' Nambu spinor, \F ) k = O^'kf) v I'j(-k)i,) T - With the 
notation introduced in Sec. [Til wc can then effectively write 



2m 



-nim, u z \ ^z.njn u z i 1 ' z .nm 



(Tz + B z 



H D = ^d n 

nm 

+ ~ ^ K [ d n (A_. nm + £?_ nm a z ) d m + d n (A^ inm + B+_ nm <r z ) d m } • 



In spinor notation, the tunneling Hamiltonian reads 



j—L^R kn 



After gauging out the phase <j) from the leads, 



T ]=L/R ^ n = cos(aL/2) e ± -^/ 4 ( % n _° 



j» 



Tj. n = =F sin(ai/2) e~ w e 



z <f>/4 ( tjn 



*j« 



After integrating out the vl/jk and taking into account the relation between the Fourier transformed 'spin-up' and 
'spin-down' Nambu spinors for the QD, see Eq. (TH)) . 

D n {w) = -j= JdTe^ T d n (T) , a x D n (-u) = -±= J dre^XW. 



we obtain the effective action 



Dm(to) 
D m (—u>) 



(Al) 



where c/^ denotes the Green's function of the uncoupled leads 



Vuj 2 + A 2 



(iw + A <7 X ) ■ 



Notice that there is no contribution from the gauged-out spin-flip terms describing spin precession along the transport 
direction, since in the second term of Eq. (|A1[) . all dd terms vanish. Finally, the action of the closed QD is 



S 



D 



T EE 



a 2 \ 



D m (w) 



+ - [D n (uj) (^4_, nm <r x + iB_ <nm <7 y ) D m (-cj) + D n (-u) (A +tnm a x - iB +<nm a y ) D m {uj) 

To simplify calculations, we now assume t* n = tj n . Introducing the multispinor ()15|) . after some algebra, we then find 

the effective action ((16)) with Eq. (fT?]) . where is the full Green's function of the QD for a = (except for the shift 
in energy), 

G-\ n , = (-iw + [e n - a 2 /(2m)} a z ) 5 nn , + ^ ± ^j"/ [-to + A cos(0/2)a x ] + {Tl ~ ^j"/ A wx(4>/2)a y (A2) 



v 7 ^ 2 + A 2 



Vw 2 + A 2 
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with the hybridization matrix given by Eq. (|18 
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FIG. 2: SO-induced oscillations of the critical current I c (in units of eA/h) in a single-level dot as a function of aL for 
b = (0, 0, 0.2r) T , r = 10A,T = 0.05A. The dot is taken in the transverse ground state n y = 0, with only one resonant level 
e Jlx = 0, for various n x . 
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FIG. 3: Zero-field critical current as a function of a in a two-level quantum dot. We take eo = 0, ui± = 20 meV, 6 = 0, 
A = 1 meV, L = 20 nm, m = 0.035m e , and T = 0.01 A. (These values are appropriate for Nb contacts.) Results are shown 
for two cases (where Ti2 = /TnTa), namely either Tn = 20A and 1^2 = (only one level couples to the leads, solid line) or 
T22 = Tn = 10A ('democratic tunneling', dashed line). 
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FIG. 4: Same as Fig. E3 but for e = 15 meV. 



